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Abstract 

Charmonia produced in initial hard parton scatterings during heavy ion collisions move with 
respect to the medium rather than flow with the medium. Lattice studies suggest that charmo- 
nium bound states at the rest are dissociated at T > 1.5 T c . We present results on momentum 
dependences of charmonium properties in a hot medium from lattice QCD Monte Carlo simula- 
tions. The dispersion relation of the screening mass and the change of correlation and spectral 
functions at various temperatures and momenta are discussed. 



1. Introduction 

In medium hadron properties can provide useful information about the hot medium formed 
during heavy ion collisions Ql, in particular heavy quarkonium can probe the deconfinement 
aspect of QGP and serve as a QGP thermometer yj. In recent years there have been great 
progress in understanding the properties of heavy quarkonium at rest in the medium from lattice 
QCD Monte Carlo simulations and heavy quark effective theories HQ]. However, quarkonia 
that are produced from initial hard parton collisions are moving with respect to medium since 
J/ip almost does not flow with the medium up to p± ~ 8 GeV/c J5[0]. Properties of moving 
heavy quarkonia have been investigated by modeling J /iff traveling in a finite temperature gas 
of mesons [7] and also in the weak coupling regime using effective theories Q8|, |9|] as well as in 
the strong coupling regime using the AdS/CFT analogy 11 Oil . These studies all suggest that the 
bound state peaks get attenuated 1F7L [8L fToll rather than unmodified Hill at large momenta. 

The moving J /iff with respect to a thermal heat bath can be treated as an external particle in 
an equilibrium background and its properties are thus suitable to be investigated on the lattice. 
From lattice QCD Monte Carlo simulations one can calculate various correlation functions, i.e. 
spatial correlation functions Gh(z, p_i_, u>„) and temporal correlation functions Gh(t, p). They are 
related with the spectral function p#(w, p, T) in the following 

G„(t, P ) = yexp H p. X )(/ H (0,0)4(r,x))= f g» C ° Sh ^ ~- p, (w , p, T), (D 



G h (z,Pj_,cl>„) = ^exp(-ip-x)(/ w (0,0)4(T,x)} = J ^" ex P^P: z )^ 



dw wpn(oj, p, T) 
7X a) 1 + cor. 



(2) 



where Jh is a local meson operator and has a form of qTnq- T# = 75 stands for the pseudo 
scalar channel (t] c ) and r# = % stands for the vector channel (J /if/) . r is the Euclidean time and 
restricted by the inverse temperature, i.e. tT < 1, z is the spatial distance and is not limited by 
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the temperature, p ± = (p x ,p y ) is the transverse momentum, p = (p_i_, a»„), <d„ = 2nnT are the 
Matsubara frequencies and x = (x,y, t). 

All information on quarkonium properties are enclosed in the spectral function and thus the 
task in principle is straightforward, i.e. to extract the spectral function from temporal or spatial 
correlation functions. Normally one extracts spectral functions from the temporal correlation 
functions as obviously the spatial correlation function has a more intricate relation with the spec- 
tral function than the temporal one. 



2. Lattice QCD results 



The charmonium properties at vanishing momentum has been recently investigated on large 
and fine quenched lattices in Ref. lll2ll . in which both P wave states (tj c and J lip) and S wave 
states (XcO and^ e i) are found to dissolve at T > 1.5 T c . In this proceeding we report results 
on the momentum dependences of charmonium properties measured on the finest lattice used in 
Ref. |H1. 

At nonzero "momentum" (p ± + or a>„ + 0), the exponential decay at large distance of the 
spatial correlator may be described by an energy E sc 



G(z, p±, (ii n ) ~ exp(-£. 5e z), E Si 



Pi 



A 2 



+ M. 



(3) 



where m sc is the screening mass which can differ from the pole mass if A(T) + 1. It is worth 
noting that the above ansatz is based on the dispersion relation in the continuum limit. 
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Figure 1 : Left: the dispersion relation of the screening mass in the pseudo scalar channel. The solid lines are linear fits 
to the data. Right: ratio of E sc at T > T c to those at T = 0.73T C with same momenta. The solid lines are results from 
free field theory with charm quark mass m c = 1.1 GeV. 



We show the dispersion relation of the screening mass in the pseudo scalar channel in the left 
plot of Fig.Q] The results are obtained from calculations performed on 128 3 x A^ T lattices at 0.73, 
1.46, 2.20 and 2.93 T c . The solid lines denote the dispersion relation obtained by fitting with 
a linear ansatz. We found that the slope of the dispersion relation almost does not change with 
temperature and only the intercept, i.e. screening mass at rest, changes with T. In the right plot 
of Fig. [1] we show the ratio of E sc atT>T c to those at T = 0.73T C with same momenta. The 
solid lines denote the results from free field theory using a charm quark mass of 1 . 1 GeV lfl2ll . At 
vanishing momentum, the deviations from unity are around 10% at the lowest temperature and 
more at higher temperature. This magnitude of deviation is comparable with those obtained from 
full lattice QCD simulations 1 14]. The ratio becomes closer to unity with increasing momentum 
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Figure 2: Ratios of longitudinal Gl(t, T, p) and transverse Gr(r, 7/, p) components of the vector correlation functions 
Gv u (t, T, p) to to the averaged correlator Ga ve(t, T, p). at 0.73, 1 .46, 2.20 and 2.93 T c . Gv u (r, T, p) are obtained using 
Jh = qy\<l and G,iv£(t, 7/, p) = Gv (t, 7", p)/3 = (Gl(t, 7\ p) + 2Gt(t, T, p))/3. Longitudinal or transverse components 
are defined such that the momentum are parallel or perpendicular to the y,, e.g. with p=(0.93,0,0) GeV, the longitudinal 
component is Gv, and the transverse components are Gy 2 and G\/, . In each plot the six sets of data points from top to 
bottom correspond to G L /G AVE at p=0.93, 1.86 and 2.97 GeV and Gj /Gave at p=2.97, 1.86 and 0.93 GeV. The solid 
lines correspond to results obtained from free field theory with charm quark mass m c =1.1 GeV. 



as in the infinite momentum limit the ratio approaches unity. It is also intuitively expected that 
the ratio gets closer to the value obtained from free field theory at higher temperatures. 

In Fig. |2] we show the ratios of the transverse and longitudinal components to the full vector 
correlator. The longitudinal components are always larger than the transverse components at a 
certain p and T. The ratios can be reproduced by the free field theory with m c = 1.1 GeV at short 
distances 111 311 . 
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Figure 3: Ratio of the measured correlator to the reconstructed correlator in the longitudinal (left) and transverse (middle) 
components in the vector channel and in the pseudo scalar channel (right) at finite momenta. 



In Fig. [3] we show the ratio of measured correlation function to the reconstructed correlation 
function, which is defined as 



;i cosh(w(r-l/2r)) / t \ 



(4) 



G rec (r,T,p) /to gg^ffl ^r^ 0.73^)' 

In practice one does not necessarily need to evaluate p(a>, T - QJ3T C , p) to get G lec (T, T, p) 
and rather can use Eq.(27) in Ref. IU2I1 to obtain G Tec (r, T,p) directly from the measured corre- 
lator G(t, T - Q.13T c ,p). One interesting thing to see is that the ratio Gl/Gl^bc almost does not 
change with momentum at each temperature while in the transverse components of the vector 
channel and in the pseudo scalar channel the momentum effects bring ratios up. This may be 
an indication of different non-trivial small frequency behavior as suggested in Ref. [15]. Note 
that calculations in Ref. lfl5ll are done for the light quark modes thus it would be more suitable 
to have a follow up study on the momentum dependences of the light correlation functions after 
Ref. iH]. 

In Fig. |4]the spectral functions extracted using the Maximum Entropy Method lfl7ll in the 
transverse components of the vector channel and in the pseudo scalar channel are shown. The 
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Figure 4: Momentum dependences of spectral functions in the vector (left) and pseudo scalar (right) channels at T = 
0.73T C and 1.46 T c . 



bands correspond to the statistical uncertainties and the horizontal error bars give the uncertain- 
ties of the peak locations at 0.73 T c . The change of the peak locations at T < T c in both channels 
are well aligned with the study on the dispersion relation of screening masses. At 1 .46 T c spectral 
functions in both channels get attenuated in the small a> region and the threshold moves to the 
large a> region. 



3. Summary 

We have studied the momentum dependences of charmonium properties by investigating on 
the temporal and spatial correlation functions. As the static r\ c and J/if/ are found to dissolve 
already at T > 1.5 T c lfl2ll . to study the properties of charmonium moving respect to the heat 
thermal bath, Monte Carlo simulations at temperature values below 1 .5 T c are crucially needed 
to examine the sequential dissociation picture of charmonia system as well as the possible mod- 
ifications of charmonium bound states at finite momentum [18]. 
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